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Abstract We develop a statistical physics theory for solid-solid phase transi- 
tions in which a metamaterial undergoes longitudinal contraction in response 
ij \ to increase in external tension. Such transitions, which are forbidden in ther- 

modynamic equilibrium, have recently been shown to be possible during the 
decay of metastable, super-strained states. We present a first-principles model 
to predict these transitions and validate it using molecular dynamics simu- 
lations. Aside from its immediate mechanical implications, our theory points 
to a wealth of analogous inverted responses, such as inverted susceptibility 
or heat-capacity transitions, allowed when considering realistic scales. 
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1 Introduction 

^ . For a closed thermodynamic system in equilibrium, the volumetric compress- 

ibility K = — yjfp is strictly non- negative because of stability constraints. 
Such constraints are manifest through first-order phase transitions and phase 
L^ ' separation via the Maxwell construction that effectively avoids any branch in 

the PV diagram that would correspond to negative k. Furthermore, during 
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Fig. 1 Inverted compressibility transitions in a 3D material in response to (a) uni- 
axial, (b) biaxial, and (c) isotropic tensional stresses (arrows). The external solid 
lines indicate the boundaries of the super-strained material prior to the increase in 
stress and the internal solid lines indicate the post-transitional state. In each case, 
the material contracts along the same directions as the increasing applied tension 
during the transition. The material and simulations are defined in the text. 



any equilibrium phase transition, the occupied volume always decreases as 
pressure increases and vice versa [T]. This analysis implicitly assumes that 
stress changes infinitely slowly. However, when the stress changes at a finite 
rate, it is common for a solid [2] (and also for polymers [3]) to occupy a 
metastable configuration for an extended period before undergoing a phase 
change, and hence the previous stability constraints no longer apply. We 
posit that, if an increase in tension causes the previous equilibrium state to 
become metastable, the solid may undergo a transition in which the decay 
of this super-strained state results in a decrease in volume as the tension 
increases. We refer to this process as an inverted compressibility transition 
(ICT). The existence of such transitions has been numerically demonstrated 
in Ref. [3] for ID and 2D systems. In this paper, we model from first princi- 
ples a related class of metamaterials that exhibit ICTs in any dimension, as 
simulated in Fig. [TJ 



2 Model System 



To develop our theory, we consider a crystal defined on a regular lattice 
in D dimensions. The constituents are shown in Fig. EJa). They consist of 
four collinear particles interacting through a potential energy V(x, y, h; a) = 
V x (x) + V y (y) +V z (y — x) + V x (h — y) + V y (h - x) - ah, where x,y, and h are 
inter-particle separations, V x , V y , and V z are pairwise interaction potentials, 
and a is the applied tensional stress. The two external particles correspond 
to lattice points, while the two internal ones correspond to edges of the 
lattice. We consider potentials such that, depending on a in the lattice, the 
constituents may occupy one of two configurations provided the temperature 
T is not too high. When a is small, a weak Vz-bond forms between the 
internal particles of each constituent, resulting in a coupled state; when a is 
increased, this bond dissociates through a decoupling transition, resulting in 
a decoupled state. This can be easily achieved with nonlinear interactions. 
We note that nonlinear interactions are receiving increasing attention in the 
context of mechanical materials — e.g., in the study of solitons [5], particularly 




Fig. 2 (a) Constituent of materials that undergo ICTs. The inter-particle distances 
and potentials are indicated on the left and right of the arrows, respectively, (b) 
Spring-string system in which the weight W is hung by two identical springs, one 
linking string, and two identical supporting strings that initially carry no tension. 
The weight may rise when the linking string is cut. 



in applications that manipulate acoustic response [HIE], and in the study of 
biomechanical networks [8j|9]. 

During a decoupling transition of the system in Fig. [5Ja), an increase in 
the applied stress causes a switch from one equilibrium (x,y,h) = (x* , y* , x* + 
y*) to a different equilibrium, (x*,y*,x* + y*). Since the derivative of the 
potential energy V{x,y,h]a) with respect to h vanishes at both equilibria, 
it follows that 

V x '{x*) + V y '(y*) = V x '(x*) + V v 'm- (1) 

If V x ' and V y ' are monotonically increasing (i.e., spring- like) we can conclude 
that x* < x* and y* > y* (or vice versa). Then, if the potentials are chosen 
such that the decrease in x overcomes the increase in y, the transition will 
indeed be an ICT. This will follow when V y hardens just after the transition 
(like a string becoming taut), and the effect can be optimized further if V x 
hardens before the transition. Specific choices for the potentials used in our 
numerics and that result in ICTs are presented in Sec. [SJ 

The decoupling transition that occurs when a is increased can be inter- 
preted more concretely in terms of the spring-string system introduced by 
Joel E. Cohen and Paul Horowitz [TU] and reproduced in Fig.^b): when the 
linking string is cut, the initially slack support strings become taut and the 
mass can, countcrintuitively, rise. This occurs because the springs go from 
a series configuration to a parallel configuration, causing the two springs 
to contract as the load becomes distributed between them. Initially each 
spring must hold the entire weight W, leading to an equilibrium height 



h = l + I — 5 



k > 



where k is the spring constant, In is the length of 



the unstretched springs, I is the length of the support strings, and S is the 
amount of slack the strings have. Once the linking string is cut, each spring 
must hold only half of the weight W, which leads to an equilibrium height 
h = In + I + |r ■ Thus the mass rises by Ah = §r — S when the linking string 
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Fig. 3 (a) A two-dimensional material consisting of a square lattice of the con- 
stituents in Fig. [U[a) for zero applied stress. Same system (b) just prior to and (c) 
just after an ICT, where the red arrows indicate the applied force profile and black 
(orange) dots indicate coupled (decoupled) states. The material touches the blue 
reference line in (b) but is entirely below it in (c), even though the applied tension 
was increased. The insets are magnified to illustrate this. When the applied tension 
is relaxed back to zero, the system returns to its initial, rectangular configuration 
(Color figure online). 



is removed, which is positive provided that S is sufficiently small. This be- 
havior has long been known [TU]. Nevertheless, likely owning to the fact that 
it involves an irreversible change of the system, this process has not been 
widely considered as a mechanism to design materials exhibiting mechanical 
behavior not found in natural materials. The ICTs undergone by the system 
in Fig. ^a) , on the other hand, only involve reversible changes and can be 
repeated by only varying the external force a, as illustrated in Fig. [3] for a 
two-dimensional lattice of such constituent systems. 



3 Equilibrium Properties 



The equilibrium properties of the proposed materials can be derived from 
the Hclmholtz free energy, A = U — TS, where U is the internal energy and 
S is the entropy To determine A given some strain e, taken as deforma- 
tion relative to zero-stress linear size, we first analyze the low temperature 
limit. For the materials we consider, ignoring any zero-point kinetic energy, 
to lowest order the free energy equals the sum of the potential energies of the 
constituents, A = Vtot- The thermodynamic equilibrium is then the state in 
which every constituent occupies the lowest-energy mechanical equilibrium 



that conforms to the external constraints. Any other mechanical equilibria 
are local minima of the free energy, corresponding to metastable states. Given 
that only the coupled equilibrium exists for small stresses and only the de- 
coupled equilibrium exists at large stresses, it follows that the coupled and 
decoupled equilibria have equal potential energy at some intermediate stress. 
Equilibrium statistical mechanics predicts that a phase transition will occur 
at this point as either the stress or the strain is varied. However, ignoring 
zero-point quantum fluctuations, the lifetimes of metastable states go to in- 
finity as temperature goes to zero. 

Further insight into this problem is obtained by considering the Debye 
model PP, which is rigorously correct in the low temperature limit and al- 
lows us to explicitly account for the effect of temperature. It is convenient to 
initially assume that all constituents are in either the coupled or decoupled 
state. The constant-strain specific heat capacity of a crystal is asymptoti- 
cally c e = aks (y) , where Ub is the Boltzmann constant, is the Debye 
temperature, proportional to the speed of sound waves in the system, and 
a is a numerical constant. Using this to determine S and U, we find that 
A = Vtot — d(d+i) (§) (^bT), which now accounts for the internal energy 
and the entropy of the system at nonzero temperature. This analysis breaks 
down only near spinodal points (i.e., points at which the coupled or decou- 
pled metastable states cease to exist), where any deviation away from T = 
will lead to a decay of the metastable states. At all other points 0^0, and 
hence the free energy is well-approximated by just the potential energy for 
sufficiently small T. This implies that the local minima of the free energy 
for small T > will correspond to those of the potential energy, so that the 
metastable states of pure phase indeed persist for small T > 0. 

For states of mixed phases, with proportion \ in the coupled state and 
proportion 1 — \ m the decoupled state, the free energy contains contributions 
from each of the separate phases plus an additional contribution due to the 
entropy of mixing, 
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v tot ,c - D (S+i) {■§:) ( k zT) 



+ (i-x) 



V tot ,d - D( S+i) \w d ) ( fcsT ) 



+ A; B T(xlogx+(l-x)log(l-x)), (2) 

where the subscripts c and d indicate the values in the coupled and decoupled 
states, respectively. The thermodynamic equilibrium is the state that mini- 
mizes A with respect to v, which for small T is y = , c , ',, ,\° / t ' d ^ \°/' c \,L\ 
to leading order. Since \ w iH now vary continuously from nearly 1 to nearly 
as the stress is increased, there is technically no equilibrium phase transition 
for any T > 0. However, if this change occurs over a scale Aa which is below 
experimental resolution, the behavior will be indistinguishable from a first 
order phase transition. Furthermore, as noted before metastability persists 
for small T > 0, so that when the applied stress or strain is varied at a finite 
rate, the response is hysteretic, just as in an ordinary first-order phase tran- 
sition. We will assume that T is sufficiently small that this is the case and 
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Fig. 4 Hysteresis loops for (a) materials exhibiting ICTs and (b) the cubic model. 
Solid lines are the thermodynamic equilibria, color dashed lines are metastable 
pure states, and black dashed lines are unstable states. The gray area includes all 
possible metastable phase mixtures, and arrows indicate transitions at spinodal 
points (horizontal when stress is controlled and vertical when strain is controlled) 
(Color figure online). 



will abuse language by referring to the transition as a "phase transition" for 
T > throughout the paper. 



4 Inverted Compressibility Transitions 



Figure 0] contrasts the hysteresis loop for materials that undergo ICTs (of 
which the materials considered here are a specific case) with the commonly 
used cubic model for solid-solid phase transitions, which does not exhibit 
ICTs. Following the figure, suppose the system begins in the lower left equi- 
librium (red solid line) as the stress is increased from small values. As the 
stress increases, it becomes possible to occupy a "flag" of metastable mixed 
states, and the system will drift into the gray area at a rate that depends 
on temperature and sample size, as quantified below. When the temperature 
is sufficiently small, the state of the system remains very near the left pure 
metastable state (red dashed line) until it reaches the vicinity of the spinodal 
point at the top of the gray area. At this point, the metastable state must 
decay, as indicated by the arrows. For the proposed materials [Fig.HJa)], this 
transition is an ICT, so that the strain decreases as the stress increases (top 
horizontal arrow), while for the cubic model [Fig. SJb)], the strain increases 
as the stress increases. When it is the strain that is controlled instead, it 
is also easy to deduce a force amplification phenomenon in which increase 
in strain leads to a transition to large restoration stress in these materials 
[top vertical arrow in Fig. SJa)], which is not present in the cubic model or 
conventional materials [top vertical arrow in Fig. HJb)]. 

We verify the above conclusions quantitatively by analyzing the decay of 
metastable states, following the general approach in solid-solid phase transi- 
tions of modeling nucleation by thermal activation [TT]. We focus on a cell 
of our material's Z?-dimensional lattice that has N constituents. For simplic- 
ity, we consider the case of forces applied equally along all D dimensions of 



the lattice. The most important macroscopic timescale is the volume fluctu- 
ation rate iy, and the temperature scale above which metastable states are 
generally short-lived is To (see Appendix). We non-dimensionalize times and 
temperatures by expressing them in units of these references. We note that 
for individual constituents, the relevant timescales are much shorter than ry . 



5 Statistical Physics Model 

5.1 Microscopic Dynamics 

To model the transitions between coupled and decoupled states, we focus 
on an individual constituent near a mechanical equilibrium (x*,y*, h*), and 
treat the remainder of the lattice as a heat bath in thermal equilibrium with 
this constituent. The constituent is assumed to have a fixed length h = h* , 
which corresponds to a canonical ensemble for its interaction with the heat 
bath. Figure [S] illustrates the vibrational modes of such a constituent in 
which the central particles move in opposite directions. These modes are of 
particular interest because they are precisely the vibrations that give rise 
to the coupling and decoupling transitions of the constituents, and thus the 
frequencies of these modes indicate how often such events may occur. 

To calculate the frequencies of these modes, we linearize the equations 
of motion assuming that the displacements from equilibrium Sxl and 8xt 
are small (Fig. [5]). We find that the respective angular frequencies of the 
longitudinal mode and transverse modes are 



^ = \ 1 , (3) 
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«r = V- " 1C ^ L , (4) 

V 777 

where 777 is the (assumed to be common) mass of the constituent particles; 
the F's are the first derivatives and the k's are the second derivatives of the 
corresponding potentials evaluated at the mechanical equilibrium point. We 
note that from the viewpoint of the entire lattice, such oscillations appear as 
optical modes of vibration. By focusing on individual constituents we have 
effectively restricted attention to optical vibrations with zero wavenumber to 
estimate the microscopic timescales. 

We assume that the timescale for which stress or strain is varied is much 
larger than Ty. Under this condition, the quasi-static approximations are 
valid and the material can be assumed to be always near (metastable) equi- 
librium, except during the transition events. In order for a constituent to ac- 
tually transition between the coupled and decoupled states, thermal fluctua- 
tions must be sufficiently large to overcome the free energy barrier separating 
the metastable state from the thermodynamic equilibrium state. We approxi- 
mate this free energy barrier as twice the potential energy differences between 
the intermediate unstable equilibrium [black dashed line in Fig.|4|a)] and the 
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Fig. 5 Modes of oscillation leading to coupling and decoupling transitions: (a) 
the longitudinal mode and (b) the D — 1 transverse modes in which the central 
particles of the constituents move in opposite directions. 



coupled (2V b c ) or decoupled (2V b d ) stable equilibrium of the constituents (as- 
suming they are initially coupled or decoupled, respectively). The factor two 
arises because the kinetic and potential energy contribute equally to the total 
energy, so only £|— of the energy is typically distributed as potential energy. 
This slightly subtle factor of two was verified as an intermediate step in our 
numerical simulations below. 



5.2 Rate of Transition Events 



We now calculate the expected rate of events that overcome the free energy 
barrier. The probability that thermal fluctuations a constituent exceed an 
energy of 2Vb is exp[-rj;], as determined by the canonical ensemble. The 
frequency at which the energy barrier is approached is derived from the 
frequencies in Eqs. ©-(HI)- During a decoupling transition, both transverse 
and longitudinal vibrations strain the V z bond, so we take the effective rate 
of energy barrier approach as [lul + {D — 1)ujt\- For a coupling transition, we 
must have a coordinated approach in both the longitudinal and transverse 
directions in order for the central particles to come together, so we use a 



reduced frequency [w L 1 + (D — l)o; 2 
for each constituent is then 



The rate of decoupling transitions 
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To proceed, suppose we increase the applied tensional stress on the N- 
constituent cell from a — at a constant rate v a = da/dt. All of the con- 
stituents in the cell are initially coupled. However, decouplings become possi- 
ble when the stress increases past some er c , where the decoupled equilibrium 



appears, and necessarily occur by some a d , where the coupled equilibrium 
ceases to exist. Then, applying our estimate for the transition rate for each 
constituent, the expected number of decoupling transitions that will occur 
by the time the stress reaches some a is given by the integral of the rate over 
time 

n d (a; N, Va , T)=N f f d (a(t>),T)dt> = N f f d (a>, T) — , (7) 

Jo Jo- v <y 

where a c < a < a d . In the second equality, we have changed the integration 
variable from time to stress and used the fact that dt' = do 1 jv a . Equation ([7]) 
should be interpreted as a valid estimate of the first transition event, which is 
all we need to complete our model for the decay of metastable states. Specifi- 
cally, we note that the transition events are largely irreversible (they tend to 
go only from metastable to stable states) and different transition events are 
not independent in the material: a single coupling or decoupling transition 
produces a cascade of further transitions, and macroscopic metastability can 
be completely and suddenly lost. Therefore, we take Eq. (|7|) as an estimate 
for the beginning of the decay of metastability. The time evolution of the 
phase transition would require the description of a kinetic relation [Till 12] . 
but such details are not necessary here. So long as n d (a; •) <C 1, we can as- 
sume that no decoupling events will occur in the cell as the system is stressed 
from a c to a. 

Taking all this together, we model the stress at which metastability is 
lost, a*, through the implicit equation 

n d (**;N,v a ,T) = l, (8) 

which we use to quantify the ICT effect. By replacing f d with f c and/or a* 
with e* and appropriately changing the limits of integration in the model 
above, similar equations can be derived for coupling transition and also for 
strain-controlled processes. 



6 Model Validation and Application 

We have verified that our model accurately predicts the onset of the decou- 
pling transitions, and ICTs in particular, using molecular dynamics (MD) 
simulations. The MD simulations were carried out using Nose-Hoover chains 
[T5] and a generalization of the Andersen barostat [T3]. Incidentally, Figs. [2 
and [3] were generated using similar MD simulations. For details on our im- 
plementation of the MD simulations, see Ref. [3]. 
In our simulations, we use the potentials 

y ( ) = J \ kx » ( x ~ - To ) 2 ' for x - Xl > (q) 

IjKMci - k Xo ) (x-xx) 2 + \k XQ (x-x f, £otx>Xi, 

v ( y ) = H k y° ( y ~ y °f > 2 2 tev^vi' (10 ) 

1 U k vi - k vo) (y - v^f + i k vo (y - vof » for y > 2/1 ■ 
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Fig. 6 Verification of our model for (a) stress- and (b) strain-controlled processes 
in a 2-dimensional square lattice. Lines represent predictions of Eq. (3) and its e*- 
analog, while points and error bars (shown when bigger than the points) represent 
averages and standard deviations of 10 MD simulations for a 20 x 20 cell. The 
different curves correspond to T = 10~ 4 (solid lines), T = 10 -3 (dashed lines), and 
T = 10 -2 (dotted lines). The strain change, At, and stress change, Aa, during 
the transitions are indicated next to the lines. Bold face Ae and Aa correspond to 
inverted responses. 



These potentials have continuous first derivatives (forces) that arc piccewise 
linear. (We have also observed ICTs in a variety of smooth potentials not 
listed here.) The V z potential is the Lennard- Jones potential, 



V z (z) = ( 



2(1 



12 



-2 



(11) 



In order to observe ICTs in our simulations, we have chosen the parameter 



values k r 



11, K 



vxo - ", -X! - 2, xo = 2, Xl = 2.2, k yo = 0.5, k yi = 20, y - 2.225, 
3/1 = 3.3, C = 0.26, and z = 0.225. 

Figure[B]comparcs the predictions of the model [Eq. ((5J) and corresponding 
equation for e*] with MD simulations for both stress- and strain-controlled 
processes. The model quite accurately predicts the general trends of the 
mctastable decay over the wide range of temperatures and stressing/straining 
rates that we simulated. In all cases, smaller stressing/straining rate v„^ and 
larger temperature T tend to lead to earlier decoupling transitions. Moreover, 
for sufficiently low temperatures (dotted and dashed lines) , these transitions 
are ICTs [Fig. EJa)] and force amplification transitions [Fig. HJb)]. 

Remarkably, as illustrated in Fig. [71 our model predicts ICTs even for very 
large, macroscopic sizes (N = 10 23 ) over very small, macroscopic stressing 
rates (y a = 10~ 10 ) for moderate temperatures (T = 10~ 3 ). In fact, this 
persistence of ICTs can be deduced by inspection from Eqs. ([7J and (jSJ). 
ICTs will occur provided metastable states persist until the stress becomes 
sufficiently large, while larger temperature T, smaller stressing rate v a , and 
larger system size TV all tend to promote earlier metastable decay. But the 
temperature dependence in Eq. (J7J is exponential, while the dependence on 
v a and TV are only polynomial. Therefore, on inverting Eq. (|5|l. it follows that 
a* depends only logarithmically on v a [Fig. |7Ja)] and N [Fig. [7Jb)], so that 
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Fig. 7 Strain change, Ae, during decoupling transitions in an isotropically stressed 
3D cubic lattice: dependence on (a) stressing rate, (b) system size, and (c) temper- 
ature. In (a) and (b), the different curves correspond to T — 10~ 5 (solid lines), 



T = 

N = 



10 (dashed lines), and T — 10 (dotted lines). When not specified, 
The shaded regions correspond to ICTs. Note that, 



10^ J and v a = 10" 



because the strain is defined relative to the zero-stress size, 100 x Ae is a percent- 
age change. 



small temperatures can easily offset their effects [Fig. Etc)]. This provides 
evidence that ICTs arc possible for a wide range of conditions. 



7 Final Remarks 



It is instructive to compare phase transitions exhibiting inverted compress- 
ibility with ordinary phase transitions, which do not. In ordinary first-order 
phase transitions, the curve of spinodal points in the phase diagram lies 
strictly within the region of phase coexistence. The emergence of an ICT im- 
plies a violation of this condition, as follows. For metastable coupled states 
reached while increasing c, the stress is larger than the stress of any state 
in the phase coexistence region (see Fig.H]). Thus, by convexity, the thermo- 
dynamic equilibrium to which these metastable coupled states decay must 
have a strain larger than the strain of any state within the phase coexis- 
tence region. Therefore, if these decays correspond to a decrease in the strain 
(as observed in ICTs), the spinodal point at which the metastable coupled 
states cease to exist must occur at a yet larger strain, and hence the curve 
of spinodal points cannot be confined to the region of phase coexistence. 
Although their potential for ICTs has not been previously recognized, the 
global behavior of spinodal curves has been of wide interest in other contexts 
for many years (see, e.g., Rcfs. |15[|T5] ). Furthermore, while ICTs have not 
been experimentally observed, we propose that with the significant advances 
made in designing both electromagnetic [17] and acoustic [18] metamaterials, 
it is possible to realize engineered materials that undergo ICTs. 

Now, one may wonder why ICTs have not been previously predicted. This 
possibility was likely not anticipated because most previous continuum [THl 
120] and lattice [21 [122] models of solid-solid phase transitions have postulated 
a free energy that is a function of a single dilation or shear strain field. Such 
models necessarily result in stress-strain relations in which the stress is a 
single- valued function of the strain, as the one shown in Fig. |4jb) . However, 
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Fig. 8 Schematics of possible hysteresis loops of an arbitrary thermodynamic 
quantity exhibiting inverted transitions. The state variable q either (a) decreases 
as the conjugate force Q increases at Q C1 or (b) increases as the conjugate force Q 
decreases at Q C2 . These transitions occur during the decay of a metastable state. 



in systems like ours, additional fields arc required to model not only the total 
length h of the constituents (corresponding to the usual strain field) , but also 
the internal x and y degrees of freedom. The many-to-many nature of the 
relationship in Fig. HJa) is the result of "projecting away" these internal 
strains. More generally, this observation indicates that any effective field 
theory that can account for ICTs will necessarily require additional fields 
(e.g., internal strains) beyond the external strain. 

Ultimately, it is the global shape of the hysteresis curve that leads to the 
inverted response, and thus analogous phenomena associated with thermody- 
namically conjugated variables other than stress and strain can also emerge. 
That is, given a thermodynamic state variable q (e.g., strain, entropy, mag- 
netization, polarization, particle number) with a conjugate thermodynamic 
force Q (e.g., stress, temperature, magnetic field, electric field, chemical po- 
tential), it generally follows from the convexity of the free energy that J* > 0. 
However, as depicted in Fig. [5J if the system occupies a metastable state, an 
increase (or decrease) in Q may force the system to decay to the thermody- 
namic equilibrium. If this results in a decrease (or increase) in the observed 
value of q, the transition will exhibit an inverted response in which ^q < 
(e.g., inverted compressibility, heat-capacity, magnetic susceptibility, electric 
susceptibility, chemical hardness transitions). We thus suggest that the neg- 
ative compressibility transitions analyzed here will have analogues in other 
thermodynamic properties, which may inspire the design of yet new classes 
of metamatcrials with unusual properties. 
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Appendix 

A. Temperature scale To 

We define the temperature scale T by setting /c_bT equal to the energy 
barrier separating the coupled and decoupled minima of the potential energy 
at the tension for which the values of these two minima are equal (i.e., the 
point at which the equilibrium phase transition would occur in the limit of 
small temperature). When T is much smaller than To, the energy barrier 
will be difficult to overcome, and metastability will persist for long periods. 
When T is much larger than To, thermal fluctuations will easily overcome 
the energy barrier, and thermal equilibrium will prevail. 



B. Timescale Ty 

For a Z?-dimcnsional cell with linear size of n primitive cells, the most im- 
portant macroscopic timescale is the inverse of the rate of the most promi- 
nent fluctuations. This corresponds to the acoustic longitudinal mode with 
wavenumber |^. For brevity, consider just a one-dimensional chain of con- 
stituents (we mention the alterations needed for the TJ-dimensional lattice at 
the end of the calculation) . Wc linearize the forces acting on the particles and 
consider small amplitude vibrations. For longitudinal modes, the equations 
of motion are 

«2„, . 

k x (vj - Uj) + k y (wj - Uj) - k y (uj - Vj-i) - k x (uj - w J _i),(12) 



dt 2 
d 2 Vj 



3 v ■ 
m ~fui~ = k v( u 3+i ~ v j)~ k x (vj -Uj) + k z (wj -Vj), (13) 

k x {u j+ i - Wj) - k y (wj - Uj) - k z ( Wj - vj), (14) 



dt 2 



where the fc's are the second derivatives of the corresponding potentials; 
Uj and ttj+i are the displacements of the exterior particles in Fig. 2(a); Vj 
and Wj are the displacements of the two interior particles. Wc insert the 
wave ansatz Uj = a u exp [i (kj — Lot)], Vj = a v exp [i (kj — out)], and Wj = 
a w exp [i (kj — Lot)] to find the dispersion relation between the wavevector k 
and the angular frequency to. The resulting secular equation is 

(— 2(k x + k y ) + mLo 2 k x + k y exp (— ik) k y + k x exp (— ik) \ 

k x + k y exp(ik) —(k x + k y + k z ) + raw 2 k z = 0. 

k y + k x exp (ik) k z —{k x + k y + k z ) + mLo 2 J 

(15) 
The acoustic longitudinal mode with wavevector k = f^ then has period 



3m (k x + ky + 2k z ) 



U 2k x ky + yk x + ky)k z 
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to lowest order in 1/n. This can be expressed as Ty = 4. / T tot n : where m to t 

is the total mass of the primitive cell and k e ff is the effective spring constant 
for the primitive cell (as found from the standard expressions for series and 
parallel configurations) . When we consider the full lattice rather than a one- 
dimensional chain, the only alteration is that k e ff is the effective spring 
constant for the lattice's primitive cell (which remains unchanged for cubic 
lattices) and mtot is the corresponding mass (1 + Z)m of the primitive cell, 
where Z is the coordination number. In the numerics, we use the zero-stress 
value of Ty as a reference timescale. The optical (longitudinal and transverse) 
modes in Eqs. ©-(HI are derived through similar analyses. 
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